Ground State Pressure and Energy Density of a Homogeneous Bose Gas in Two 

Dimensions 
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We consider an interacting homogeneous Bose gas at zero temperature in two spatial di- 
mensions. The properties of the system can be calculated as an expansion in powers of g, where 
g is the coupling constant. We calculate the ground state pressure and the ground state energy 
density to second order in the quantum loop expansion. The renormalization group is used to sum 
up leading and subleading logarithms from all orders in perturbation theory. In the dilute limit, 
the renormalization group improved pressure and energy density are expansions in powers of the 
T-matrix. 



I. INTRODUCTION 



The remarkable achievement of Bose-Einstein con- 
densation (BEC) of alkali atoms in harmonic traps 
has created an enormous interest in the properties of di- 
lute Bose gases (For a review, see e.g. Ref. [Q and refer- 
ences therein). 

The homogeneous Bose gas in three dimensions has 
been studied in great detail over the past 50 years (For a 
review, see e.g. |5j). At zero temperature, the quantum 
loop expansion is essentially an expansion in powers of 
yj pa 3 , where a is the two-body S-wave scattering length 
and p is the density. Lee and Yang Q were the first to 
calculate the leading quantum correction to the energy 
density. Part of the second quantum correction to the 
energy density was obtained by Wu, by Hugenholz and 
Pines, and by Sawada Q . Only recently a complete two- 
loop result has been obtained by Braaten and Nieto ||. 
The result depends, in addition to the scattering length, 
also on the scattering amplitude for 3 — > 3 scattering. 

The homogeneous Bose gas in two dimensions has 
also been studied extensively. The chemical potential 
and ground state energy density of a two-dimensional ho- 
mogeneous Bose gas were first calculated by Schick [ [To|| . 
By summing up the ladder diagrams contributing to the 
chemical potential, he was able to show that, in the di- 
lute limit, the relevant expansion parameter is not the 
coupling constant g, but rather the effective interaction 
[log(pa 2 )] _1 , where p is the density and a is the range 
of the interaction. The expansion parameter is the two- 
body T-matrix, and Schick determined the leading-order 
results in this expansion. Corrections to these results 
have been considered by several authors pd[-fiq|. Very 
recently, a formal proof of the result by Schick was given 
by Lieb and Yngvason jl7|,[l8) ■ 

In the present paper, we reconsider the homogeneous 
Bose gas in two dimensions at zero temperature. We 
calculate the pressure and energy density of the ground 
state to second order in the quantum loop expansion. We 



also apply the renormalization group to sum up leading 
and subleading logarithms from all orders of perturbation 
theory. In the dilute limit, the renormalization group 
improved pressure and energy density are essentially ex- 
pansions in powers of the T-matrix. 

The outline of the paper is as follows. In Sec. II, we 
briefly discuss the perturbative framework developed in 
Ref. H to calculate the ground state properties of a ho- 
mogeneous Bose gas. In Sec. Ill, we calculate the ground 
state pressure to two loop order. In Sec. IV, we calcu- 
late the ground state energy density through two loops. 
In Sec. V, we apply the renormalization group to sum 
up leading and subleading logarithms from all orders in 
perturbation theory. Finally, we summarize our results 
in Sec. VI. Calculational details are included in an ap- 
pendix. 



II. PERTURBATIVE FRAMEWORK 

In this section, we discuss the perturbative frame- 
work set up in Ref. [pi to calculate the effects on the 
ground state from quantum fluctuations around the mean 
field. 

The action is 




d 2 x V*(x, t) 



V*(x,i) 



d z x \ d z x' 



^(x,i)V(x',f) 



7/>*(x,i)V*(x',i)V (x-x') x 



(1) 



i/j*(x,t) is a complex field operator that creates a bo- 
son at the position x, /i is the chemical potential, and 
Vb(x) is the two-body potential. In the following, we set 
h = 2m = 1. Factors of h and 2m can be reinserted using 
dimensional analysis. 
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The action Eq. ([!]) is symmetric under a phase trans- 
formation 



V>(x,t) -> e iQ ^(x,t) 



(2) 



The [/(l)-symmetry ensures that the density p and 
current density j satisfy the continuity equation 



V-j = 



(3) 



In the ground state, the current density j vanishes iden- 
tically and the condensate has a constant phase. The 
[/(l)-symmetry can then be used to make the conden- 
sate real everywhere. 

If the interatomic potential Vo(x) is short range, it 
can be mimicked by local interactions. If the energies 
are low enough, the scattering amplitude can be approx- 
imated by s-wave scattering and the action Eq. (^) can 
be approximated by fcjj 



S 



dt / d 2 x tp* 



,d_ 

'dt 



(4) 



Here, g is a coupling constant that must be tuned to re- 
produce some low-energy observable of the true potential 
Vb(x). 

The quantum field theory defined by the action 
Eq. (|4|) has ultraviolet divergences that must be removed 
by renormalization of p and g. There is also an ultravi- 
olet divergence in the expression for the density p. This 
divergence can be removed by adding a counterterm dp. 
Alternatively, one can eliminate the divergences associ- 
ated with p and p by a normal-ordering prescription of 
the fields in Eq. (^). The coupling constant is renormal- 
ized in the usual way by replacing the bare coupling with 
the physical one. 

If we use a simple momentum cutoff M to cut off 
the ultraviolet divergences in the loop integrals, there 
will be terms proportional to M p , where p is a positive 
integer. There are also terms that are proportional to 
logM. The coefficients of the power divergences depend 
on the regularization method and are therefore artifacts 
of the regulator. On the other hand, the coefficients of 
log(M) are independent of the regulator and they there- 
fore represent real physics. In this paper, we use dimen- 
sional regularization to regulate both infrared and ul- 
traviolet divergences. In dimensional regularization, one 
calculates the loop integrals in d = 2 — 2e dimensions 
for values of e where the integrals converge. One then 
analytically continues back to d = 2 dimensions. With 
dimensional regularization, an arbitrary renormalization 
scale M is introduced. This scale can be identified with 
the simple momentum cutoff mentioned above. An ad- 
vantage of dimensional regularization is that it automat- 
ically sets power divergences to zero, while logarithmic 
divergences show up as poles in e. In two dimensions, 



the one-loop counterterms for the chemical potential p 
and the density p are quadratic ultraviolet divergences, 
while the one-loop counterterm for the coupling constant 
g is a logarithmic ultraviolet divergence. At the two-loop 
level, the counterterms for the chemical potential and the 
density are also quadratic divergences. The counterterm 
for the coupling constant is a double logarithmic diver- 
gence. 

We next parameterize the quantum field ip in terms 
of a time-independent condensate v and a quantum fluc- 
tuating field ip: 



tp = v + ip . 



(5) 



The fluctuating field can be written in terms of two real 
fields: 



i> = 4| + #2) 



(6) 



Substituting Eq. (^) into Eq. the action can be de- 
composed into three terms 



S[lp] = S[v] + S'freelV'l!^] + S int [v , ^ , 1p 2 ] 

S\v] is the classical action 



S[v] 



dt / d 2 x 



1 



fiv 



(7) 



(8) 



while the free part of the action is 



dt / d 2 x 



+\ipi (V 2 + X) tpi + X -i> 2 (V 2 + Y) ij 2 



The interaction part of the action is 



»S'int[f,'0i ; ' ! /'2] = I dt I d 2 x V2T^i 



(9) 



The sources in Eq. (10) are 



T = [a - gv 2 ] v 

X=[n- 3gv 2 ] 

Y=[p- gv 2 ] , 
Z = -gv . 



(10) 



(11) 
(12) 

(13) 
(14) 



The propagator that corresponds to the free action 
Sfteeliputpz] in Eq. @ is 



p 2 -Y 



Jul 



p 2 -X 



(15) 
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Here p is the wavevector, u> is the frequency, and e(p) is 
the dispersion relation: 



e(p) = ^/(p 2 -X)(jp 2 -Y) 



(16) 



The value of the condensate vq that minimizes the clas- 
sical action is given by the equation T = 0. Both the 
propagator Eq. ( |l5| ) and the dispersion relation Eq. fll6| ) 
greatly simplify at t> , since Y — here. The dispersion 
relation then becomes gapless, which reflects the spon- 
taneous breakdown of the £/(l)-symmetry (Goldstone's 
theorem). The dispersion relation is linear for small 
wavevectors and is quadratic for large wavevectors, which 
is that of a free nonrelativistic particle. The propagator 
is defined with an ie prescription in the usual way. The 
diagonal parts of the propagator are denoted by a solid 
and a dashed line, respectively. The off-diagonal parts 
are denoted by a solid-dashed and dashed-solid line, re- 
spectively. 

The partition function Z can be expressed as a path 
integral 



(17) 



All the thermodynamic observables can be derived from 
Z. For instance, the free energy density T is given by 



; log£ 
' VT ' 



(18) 



where VT is the spacetime volume of the system. 

The density p is given by the expectation value 
(i/)'ip) in the ground state. It can therefore be expressed 
as 



pip) = - 



dp 



(19) 



The energy density S is given by the Legendre transform 
of the the free energy density 



£(p) = J~(p) + PP 



(20) 



it can be shown that all the one-particle reducible graphs 
(those that can be disconnected by cutting a single line) 
vanish. We then have 



F(p) = Qo(p,v) + fii(/x,u) + Q 2 (p,v) + 
The condition Eq. ( |22] ) is equivalent to 

£-0. 

ov 



(23) 



(24) 



The loop expansion Eq. (|23|) does not coincide with the 
expansion in powers quantum corrections. To obtain the 
expansion in powers of quantum correction, we must ex- 
pand the condensate v about the classical minimum vq: 

v = v + Vi + v 2 + ••• • (25) 

By substituting Eq. ( [2l|) into Eq. (pif ) , we obtain the first 
quantum correction v± to the classical minimum: 



dQi(p,v) 



dv 



d 2 n a (p,v) 



dv 2 



(26) 



v=v l 

This first quantum correction to the free energy density 
is 



Fi{p) = HiCMjVo) , 



(27) 



and the second quantum correction to the free energy 
density is 



T 2 ip) = Q, 2 ip,v ) + vi 



d£liip,v) 



i 2 d 2 n ip,v) 



dv 2 



dv 



v=v 



(28) 



III. PRESSURE TO TWO LOOPS 

In this section, we calculate the pressure as a func- 
tion of the chemical potential p and the renormalized 
coupling g to two loops. 



At this point it is convenient to introduce the ther- 
modynamic potential Qip,v). The thermodynamic po- 
tential is given by all one-particle irreducible vacuum 
graphs and can be expanded in the number of loops 



Cl(fj,,v) = Q iP,v) + Q,iip,v) + U 2 {p,v) 



(21) 



where the subscript n indicates the contribution from 
the nth order in the loop expansion. The free energy T 
is given by all connected vacuum graphs and is indepen- 
dent of the condensate v. If we evaluate at a value of 
the condensate that satisfies the condition 



(22) 



A. Mean-field free energy 

The thermodynamic potential in the mean-field ap- 
proximation is 



ft (fj,,v) = -pv 2 + ^gv 4 



(29) 



The mean-field free energy is given by the classical ther- 
modynamic potential evaluated at the classical minimum 
v : 

,2 



Faip) 



P_ 

2ff 



(30) 
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B. One-loop free energy 



C. Two-loop free energy 



The one-loop contribution to the free energy is 
^i(At) =^ 7 i q (m) + A 1 ^, 



where 

„ , . i f du f d d p , , „ i 

^ = 2/ W W7 logdctD 

and A1J 7 is the one-loop counterterm 



(31) 



(32) 





A!^ = -^A lA1 + ^A lff 
9 2 5 2 



(33) 



The propagator D(ui,p) is evaluated at the classical min- 
imum, where Y = 0. By integrating over uj, Eq. ( p{2] ) 
becomes 



(34) 



;7 0) -i(2/j) , 



where 7 m ,„ is defined in the appendix. Using Eq. ( |A.9| ) 
in the appendix, we obtain 



167T 



i-^-5 + o(e) 



(35) 



where 7 = log (/i/2Af 2 ). The counterterms at one-loop 
are fijl: 



Ai/i = 0, 
A 5 2 



(36) 
(37) 



Adding Eqs. ( p3[ ) and (35), we obtain the one- loop con- 
tribution to the free energy 



16tt 



1 



(38) 



Note that T\ (/x) is independent of the coupling constant 
g. Adding Eqs. ( |30| ) and (|39|), we obtain the one-loop 
approximation to the free energy: 



■^0+1 (M) = 1" T^- 

(7 1D7T 



1 

L+ 2 



(39) 



where the coupling constant is evaluated at the scale M, 
9 = 9(M). 



FIG. 1. Two-loop vacuum graphs contributing to the free 
energy T. 

The two-loop contribution to the thermodynamic 
potential f2 2 (/z, i>) isgiven by the one-particle irreducible 
graphs shown in Fig. [1} The expression for the diagrams 
is 



f2 2 (/x,u) = ^7 



9_ 

32 



[3/2 x + 27_i,_ 1 Ii,i 



+3^i,-i] , (40) 

where 

J = 6Jo : o,i ^ '7-1,-1,1 — 3 Ji,i,i — 2J_i,o,o > (41) 

and integrals Ji, m , n are defined in the appendix. The 
first correction v\ is given by Eq. (Um: 



gyp 

8/x 



[3/i 



(42) 



The two-loop contribution to the free energy is then 



^2 00 



.9/' 



J 
A 2 ^ 



16 



-1 - 2/_ 1) _ 1 7i 



37 



li] 



(43) 



The counterterm A 2 .7" is 

A 2 :F = -7i,iAi/i - -A 2 /i 
2 3 



1 



A 4 



— (Aiii) z + ^AiitA^ - ^A l5 



(44) 



2g * ^' ' ^ " 2 5 3 
At the two-loop level, the counterterms are given by 

A 2A i = , (45) 
9 3 



A 2ff = 



647r 2 e 2 



(46) 



The integrals J;, m ,„ are ultraviolet divergent. After sub- 
tracting the divergent part, the remainder must be cal- 
culated numerically. We evaluate the expression in the 
appendix. The result is 



J 



-27 0)1 [7i 



7_i, 



Jn 



(47) 



4 



where J num = —3.52 x 10 5 /i. The final result for the 
two-loop contribution to the free energy is 

,2 



gp 
'64tt ; 



[1 + C] 



(48) 



Here C = 2.78 x 10~ 3 . Adding Eqs. <Mj, (B9h, and (48) 



we obtain our final result for the free energy to second 
order in the quantum loop expansion 



0+1+2 



(/') 



25 



JJ_ 
16tt 



L 



gp 

64tt 2 



[1 + C] , (49) 



D. Pressure to two loops 

The pressure V is given by —T. The pres- 
sure through two loops is given by minus the sum of 
Eqs. ©, © and @: 

,2 

[1 



9 9 

^0+1+2^)-^-^ 



1 

i+ 2 



64tt 2 



The coupling constant <? in Eq. (||) satisfies 



M q 

dM y 



C] . (50) 
(51) 

where the f3- function is a polynomial in g. Normally, 
the /3-functions are known only up to a certain order in 
the quantum loop expansion. In the present case, the 
one-loop result for the /3-function is exact and (3(g) = 
g 2 /Air [^9). From Eq. j5l|), one can easily check that our 
result Eq. ([30]) for the two-loop pressure is independent 
of the scale M up to correction of order g 2 . 



IV. ENERGY DENSITY TO TWO LOOPS 

In this section, we derive the energy density £ as a 
function of the density p and the renormalized coupling 
g to two loops. 

A. Mean field energy density 

Using Eqs. (|l^) and (^) the density in the mean- 
field approximation is 

poGO = - ■ (52) 
9 

The chemical potential is obtained by inverting Eq. (52): 

Po(p)=PP- (53) 

Using Eqs. (|2^) and (|52|), the energy density in the mean- 
field approximation is given by 

1 



£o(p) = 77.9 p 2 



(54) 



B. One-loop energy density 



Using Eqs. ( |19[ ) and (|39j) we obtain the density in 
the one-loop approximation 



/ N, P P 

P0+1(M) = ___ 



L + l 



(55) 



Inverting Eq. ( pq ) to obtain /1 as a function of p, one 
finds 



r \ ^9 2 P 
Po+i(P) = 9P+ -5— 



L + l 



(56) 



where L = log(gp/2M 2 ) and g — g(M). Using 
Eqs. (|20|), ( |39"| ) and (|56|), the energy density in the one- 
loop approximation becomes 



c ( \ 1 2 , 9 2 P 2 



L 



(57) 



This agrees with the result obtained by Lozano |2(J , and 
by Haugset and Ravndal pjj. 



C. Two-loop energy density 

Using Eqs. ( |l9| ) and (|9|), we obtain the density in 
the two-loop approximation: 



/i p 

P0+1+2W = r- 

5 8?r 



L + l 



+ #2 [1 + C] ■ (58) 



Inverting Eq. (pq ) , we obtain the two-loop expression for 
the chemical potential: 



Mo+i+2(p) = gp 



g 2 p 

8tt 



L+l 



64+ 



[L 2 +3L- 2C] 



(59) 



The final result for the two-loop energy density then be- 
comes 



£0+1+2 (p) = \gp 2 + -^-p 2 

2 1d7t 
o 3 

■ y 2 r f 2 



128+ 



p 2 [L 2 + 2L - 1 - 2C] 



(60) 



Using the renormalization group equation ( |5l| ) for the 
running coupling constant, we see that the two- loop re- 
sults for the density, chemical potential, and energy den- 
sity are independent of the renormalization scale M up 
to corrections of order g 4 . 
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V. RENORMALIZATION GROUP 

In Sees. Ill and IV, we showed that the results for the 
pressure, density, chemical potential, and energy density 
are independent of the arbitrary renormalization scale M 
that is introduced with dimensional regularization. 

The two-loop results for the pressure and density 
include a logarithm L. The renormalization group can 
be used to absorb this logarithm. The solution to the 
renormalization group equation (Ell) is 



g(M ) 



9(M) 



1 



8 



(61) 



By substituting the running coupling constant into the 
expressions for the pressure and density in the mean-field 
approximation and choosing the renormalization scale 
Mo = yj A*/2, we absorb this logarithm. In a similar man- 
ner, we can sum logs of the form g n+1 L n (n = 1, 2, 3...) 
by using the running coupling in the remaining terms in 
expressions for the pressure and density. These terms are 
generated by expanding out Eq. (^Tj) in powers of L. The 
result is 



PW = — - 5- 



90P r. 



C] , 



V{H) = 



2.9o 



A* 

32tt 



64tt 2 



[1 + C] , 



(62) 
(63) 



where go = fir (Mo). In the same way, we can sum leading 
(of the form g n+1 L n , where n — 1,2,3...) and sublead- 
ing logs (of the form g n+m L n , where n = 1,2,3..., and 
to = 2, 3...,) from all orders of perturbation theory in the 
chemical potential and energy density by using a running 
coupling constant g(Mo) and choosing the renormaliza- 
tion scale Mo = J p/2. These terms are generated by 
expanding out Eq. (|6l[). The result is 



2 

M(p) = .9oP+^[log(. 9 o) + l] 



9oP 
64tt 2 



[log 2 ( 50 ) + 31og( ff0 )-2C] , 



(64) 



log(ffo) + 2 



2 16-7T 

+ 1 §§£a [WGto) + 21og( 50 ) -1-2C], (65) 

where g = g(M ). 

We have approximated the two-body potential by 
local interactions. In real systems, however, the interac- 
tions have a finite range a and so 1 /a provides a natural 
ultraviolet cutoff M. In the dilute gas limit, we have 
|log(/ia 2 /2)| > 1 and |log(pa 2 /2)| > 1 and 
the running coupling constant is approximately equal to 
either -87r[log(/ia 2 /2)]- 1 or -87r[log(pa 2 /2)]- 1 . This is 



essentially the two-body T-matrix and so this is the rele- 
vant expansion parameter. By substituting these expres- 
sions for the coupling constant into to the expressions 
for the pressure, density, chemical potential, and energy 
density, we obtain the renormalization group improved 
results 



pip) 

Kp) 
^p) 



1 



1 + -L- 1 + 2{1 + C)L~ 2 



[1 + L- 1 + 2(1 + C)L- 2 ] , 
[1 + (log[-Z/87r] - 1) L- 1 



(66) 
(67) 



fi 2 L 

~ 167T 

fiL 

8np 
77 

+ log 2 (log[-Z/87r]) L~ 2 - 3 (log[-Z/87r]) L~ 2 
-2CL- 2 ] , (68) 



47rp 2 

z 



l + i(21og[-Z/8vr]-l)Z- 1 



+ log 2 (log[-i/87r]) L- 2 - 2 (log[-i/87r]) L~ 2 
-(1 + 2C)Z- 2 ] , (69) 

where where L = [log(/xa 2 /2)] . and Z = [log( j oa 2 /2)] . 
Eqs. (|66])-(|69|) are the main result of the present pa- 
per. The leading order results were first obtained by 
Schick |[(|, while the leading corrections have been cal- 
culated in Refs. ITll-Ol. The last terms are new result. 



VI. SUMMARY 

In the present paper, we have studied a two- 
dimensional interacting homogeneous Bose gas at zero 
temperature. We have calculated the ground state pres- 
sure and energy density to second order in the quantum 
loop expansion. The results are independent of the arbi- 
trary renormalization scale M. 

We have applied the renormalization group to 
sum up leading and subleading logarithms from all 
orders in perturbation theory. The renormalization 
group improved pressure and density are expansions in 
[log(/^a 2 /2)] , while the chemical potential and energy 

density are expansions in [log(pa 2 /2)] . These are es- 
sentially expansions in the two-body T-matrix, and we 
have obtained the ground state pressure, density, chem- 
ical potential, and energy density to next-to-next-to- 
lcading order. 
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APPENDIX A: FORMULAS 

The loop integrals that appear in our calculations 
involve integrations over the energy u> and the spatial 
momentum p. The energy integrals are evaluated using 
contour integration. 

The specific one-loop integral needed is 



du) 



1 



2tt [lo 2 - e 2 {p) +ie] 2e(p) ' 
The specific two-loop integrals needed are 
duj\ f du>2 1 



(A.l) 



2tt J 2tt [u\ -e 2 (p) + ie] 
1 

[iv 2 - € 2 {q) + i^Jl + lo 2 ) 2 - € 2 {r) + ie] 
1 

4e(p)e(g)e(r)(e(p) + e(g) + e(r)) ' 

du>l f du>2 ^\0J2 



(A.2) 



2tt J 2?r [w\ -e 2 (p) + ie] 
1 

M2 ~ e 2 {q)+ie\[{LOi+uJ 2 ) 2 - e 2 (r)+ie] 

1 



4e(r)[e(p) + e(q) + e(r)] 



(A.3) 



Here r = |p + q|. 

Some of the one-loop momentum integrals are in- 
frared divergent or ultraviolet divergent or both. They 
can be written in terms of the integral I m ,n, which is 
defined by 



Im,n 



TM 2 \ e f d d P 



-in 



P 



(2ir) d p n (p 2 + A 2 ) n / 2 ■ 



(A.4) 



Here, M is a renormalization scale that ensures that Im.n 
has the canonical dimension also for d 7^ 2. 7 w 0.5772 is 
the Euler-Mascheroni constant. With dimensional regu- 
larization, 7 m)7i is given by the formula 



z 1 M 2 
(2v) d V 4?r 

i-T 

2 



2m-2ii 



A d + 



r(fca+m)r(n-m-g) 



2r(f) 



(A.5) 



where fid — 2ir d / 2 /T[d/2] is the area of the ci-dimensional 
sphere. 

The integrals I m .n satisfy the relations 



The first relation follows directly from the definition of 
Im.n- The second relation follows from integration by 
parts, while the last is simply an algebraic relation. 

In two dimensions, these integrals have logarithmic 
and power ultraviolet divergences. The power diver- 
gences are set to zero in dimensional regularization, while 
the logarithmic divergences appear as poles in e. The 
specific integrals are 



*0,-l - 



1-1,-1 = 



A 4 
32tt 
1 

+ 2 
A 2 



L 



I- L-l 

e 2 

+ L -I 



7T 



1 



— i--L 
8tt 1 e 



(L-l) 



h,i = 



1 

+ 2 

A 2 , . 

<^ - -L-l 

8tt I e 



1 + 



l) 2 -4 



7T 

T 



e + (D(e 2 



O 



1 J 1 r 1 

47U- L+ 2 



r 2 


e + o( e 2 )} 


+ T 





(A.9) 

(A.10) 

(A.ll) 
(A.12) 



where L = log(A 2 /4M 2 ). 

The two-loop integrals needed can be expressed in 
terms of Ji m ,n 



4-7T 



d d p f d d q 



(2ir) d J (2ir) d 



p/y/p 2 + A 2 q/y/q 2 + A 2 r/Vr 2 + A 2 



Pa/p 2 + A 2 + (? V 1 ? 2 + A 2 + tVt* 2 + A 2 



(A.13) 



In two dimensions, these integrals have quadratic and 
double logarithmic divergences that cancel in the partic- 
ular combination in Eq. ([43|), leaving us with a logarith- 
mically divergent integral. We write the integral as 



J — Jdi- 



(A.14) 



where the ultraviolet divergence of the integral J has 
been isolated: 



Jdiv 2 



47T 



yV + A 2 
P 



d d p 
(2n) d 

d d q 



P 



A 2 



(27r) d q^TT 2 



(A.15) 



rfA 2 



I'm . ri ^ I; 



m+l,n+2 j 



(d + 2m- n) I m<n = nI m+2 m+2 , 

A Im.n Im — 2 Im+l.n 



(A. 6) The first term inside the square brackets in Eq. (A.15) 
has only a power divergence, and so it is set to zero in di- 

(A.7) 

mensional regularization. In terms of I m .m the remaining 
(A. 8) terms can be written as 



7 



Jd 



-2/0,1 [h 



■1-1,-1 



(A.16) [2 



The remaining finite part of the integral J can be evalu- 
ated directly in two dimensions and reads 




vV + A 2 P 



P 



A 2 



y/ p 2 + a 2 ^ 2 + A 2 Vr 2 + A 2 



r \/p 2 + A 2 \/<? 2 + A 2 
pqVr 2 + A 2 



P\/p 2 + A 2 + g^g 2 + A 2 + r\Jr 2 + A 2 



g^g 2 + A 2 

y/p 2 + A 2 
P 



2 - 



vV + A 2 



(A.17) 



Since the only scale in the integrand in Eq. (A.17) is A, 
it follows from dimensional analysis that J num is propor- 
tional to A 2 . The numerical value is 



Jn 



-1.76 x 10 A 



5 A 2 



(A.18) 
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